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Abstract. Radial velocity studies of X-ray binaries provide the most solid evidence for the 
existence of stellar-mass black holes. We currently have 20 confirmed cases, with dynamical 
masses in excess of 3 Mq. Accurate masses have been obtained for a subset of systems which 
gives us a hint at the mass spectrum of the black hole population. This review summarizes the 
history of black hole discoveries and presents the latest results in the field. 
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1. Introduction 



Theoretical black holes (BHs) have been considered for almost a century, although 
their first observational evidence was not found until quite recently, during the past 3 
decades. Since then BHs have become essential ingredients in the construction of modern 
1 astrophysics on all scales, from binaries to galaxies and AGNs. It is stellar-mass BHs 
that offer us the best opportunity to study these objects in detail. Their implications 
are wide ranging, from late evolution of massive stars, supernovae models, production 
of high-energy radiation, relativistic outflows, chemical enrichment in the Galaxy, etc. 
But still the most compelling evidence for the existence of stellar-mass BHs relies on 
dynamical arguments, which will be the focus of this review. Section 2 summarizes the 
first observations of historical BH "candidates", whereas Sect. 3 reviews the most recent 
discoveries. In Sect. 4 we deal with BH demography and their mass distribution. Finally, 
we present our conclusions in Section 5. 
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^ | 2. Early discoveries: the first BH "candidates" 

In the late 1960s X-ray detectors onboard satellites revolutionized astronomy with 
• ■ *h . the discovery of an unexpected population of luminous X-ray sources in the Galaxy. The 
energetics (with L x ~ ^Edd) together with their short timescale variability (down to 
milliseconds) lend support to an interacting binary model where X-rays are supplied by 
accretion onto a collapsed object (Shklovskii 1967J) . It is now well stablished that there are 
two main populations of X-ray binaries: the high mass X-ray binaries (HMXBs), contain- 
ing O-B supergiant donor stars, and the low mass X-ray binaries (LMXBs), with typically 
short orbital periods and K-M donors. The optical flux in LMXBs is triggered by repro- 
cessing of the X-rays into the accretion disc whereas in HMXBs it is dominated by the hot 
supergiant star (see reviews in Charles & Coe 2006 and McClin Fock fc Remillard 2006[) . 

Therefore, it was not surprising that one of the first optical counterparts to be 
identified was the 9th magnitude supergiant star HD 226868, associated with the HMXB 
Cyg X-l. But, remarkably, it showed radial velocity variations which made it a prime 
candidate for a stellar-mass BH (Webster & Murdin 1972. Bolt on 1972ft . The supergiant 
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star was shown to move with a velocity amplitude of ^64 km s~ 1 (later refined to 75 km 
s _1 ) in a 5.6 day orbit due to the gravitational influence of an unseen companion (see Fig. 
1). The orbital period P rb and the radial velocity amplitude K, combine in the mass 
function equation /(M x ) = K 3 P OTh /2TrG = AfJ sin 3 if (M x + M c f which relates the 
mass of the compact object M x with that of the companion star M c and the inclination 
angle i. The mass function f(M x ) is a lower limit to M x and, for the case of Cyg X-l 
is 0.25 M Q . The key factor here is M c which, for a HMXB, is a large number and has a 
wide range of uncertainty. If the optical star were a normal 09.7Iab its mass would be ~ 
33 M Q which, for an edge-on orbit (i = 90°), would imply a compact object of ~ 7 M Q . 
However, the optical star is likely to be undermassive for its spectral type as a result of 
mass transfer and binary evolution, as has been shown to be the case in several neutron 
star binaries (e.g.[Rappapor t fc Joss 1983 ). In fact, it could be undermassive by as much 
as a factor of 3 given the uncertainty in distance, log g and T c g. A plausible lower limit 
of 10 M Q , combined with an upper limit to the inclination of 60°, based on the absence 
of X-ray eclipses, leads to a compact object of >4 M Q (jBolton 1975|) . 
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Figure 1. Radial velocity curve of 
HD 226868, the 09.7Iab companion 
star in the HMXB Cyg X-l, folded 
on the 5.6 day orbital period. After 
Webster & Murdin (1972). 



The importance of this result rests on the fact that there is a maximum mass for neu- 
tron stars (NS) to be stable against gravitational collapse fOppcnhci mer fc Volkoff 1939[ ). 
This maximum depends on the equation of state (EoS), which is uncertain in the high 
density regime because of the poorly constrained many-body interactions. However, 



Rhoades & Rufflni (1974) showed that an upper limit of ~3.2 Mq can be derived assum- 



ing that causality holds beyond densities where the EoS starts to be uncertain which, at 
the time, was 1.7 times the density of nuclear matter p nm = 2.7 x 10 14 g cm -3 . More 
recently, a new limit of 2.9 Mq was obtained using modern EoS which are accurate up 
to 2 x p nm (Kalogera & Baym 1996). This can be further boosted by up to ~25% if the 
NS rotates close to break-up (Friedman & Ipser 1987). Therefore, the compact object in 
Cyg X-l, with M x 4 Mq, is a very strong BH candidate. 

In 1975 the satellite Ariel V detected A0620-00, a new X-ray source which displayed 
an increase in X-ray flux from non-detection to a record of ~50 Crab. It belongs to 
the class of X-ray transients (XRTs, also called X-ray Novae), a subclass of LMXBs 
which undergo dramatic episodes of enhanced mass-transfer or "outbursts" triggered by 
viscous-thermal instabilities in the disc (e.g. King 1999). During outburst, the companion 
remains undetected because it is totally overwhelmed by the intense optical light from 
the X-ray heated disc. However, the X-rays switch off after a few months of activity, 
the reprocessed flux drops several magnitudes into quiescence and the companion star 
becomes the dominant source of optical light. This offers a very special opportunity 
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to perform radial velocity studies of the cool companion and unveil the nature of the 
compact star. The first detection of the companion in A0620-00 revealed a mid-K star 
moving in a 7.8 hr period with velocity amplitude of 457 km s _1 . The implied mass 
function was 3.2 ± 0.2 M©, the largest ever measured (McClin tock fc Remillard 1986|) . 
An absolute 3<r lower limit to the mass of the compact star of 3.2 M Q was established 
by assuming a very conservative low-mass companion of 0.25 M a and i < 85°, based on 
the lack of X-ray eclipses. This exceeds the maximum mass allowed for a stable NS and 
hence it also became a very compelling case for a BH. 



3. From "candidates" to confirmed BHs 

In the 1980s there was a hot debate about the real existence of BHs. On one hand 
there were 3 strong candidates, the HMXBs Cyg X-l and LMC X-3 and the transient 
LMXB A0620-00, all with lower limits to M x very close to the maximum mass for NS 
stability. On the other, alternative scenarios were proposed to avoid the need for BHs such 
as multiple star systems (|Bahcall et al 1974)) or non-standard models invoking exotic EoS 
for condensed matter. An example of the latter are Q-stars, where neutrons and protons 
are confined by the strong force rather than gravity and can be stable for larger masses 
(jBahcall et al 1 990). In this context, it was proposed that the "holy grail in the search for 
BHs is a system with a mass function that is plainly 5 M© or greater" ()McClintock 1986]) . 




Figure 2. Radial velocity curve 
of the K0 companion in the tran- 
sient LMXB V404 Cyg during 
quiescence. This graph contains 
velocity points obtained between 



1991 and 2005. 
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In 1989, the X-ray satellite Ginga discovered a new XRT in outburst named GS 
2023+338 (=V404 Cyg). Its X-ray properties drew considerable attention because of the 
exhibiti on of a possible luminosity saturation at L x ~ 10 39 erg s _1 and dramatic vari- 
ability (Zycki, Done & Smith 1999). Spectroscopic analysis during quiescence revealed a 
K0 star moving with a velocity amplitude of 211 km s _1 in a 6.5 day orbit (see Fig. 2). 
The mass function implied by these numbers is 6.3 ± 0.3 Mq, and hence the compact 
object must be more massive than 6 M©, independent of any assumption on M c and i 
(jCasares et al 1992|) . This remarkable result established V404 Cyg as the "holy grail" BH 
for almost a decade. Since then, many other BHs have been unveiled through dynamical 
studies of XRTs in quiescence, seven others with mass functions also in excess of 5 M© . 
This has been possible thanks to the improvement in spectrograph performance on a new 
generation of 10-m class telescopes over the last decade. 

Table 1 presents the current list of 20 confirmed BHs based on dynamical arguments, 
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Table 1. Confirmed black holes and mass determinations 



oystem 


P 

J orb 


f(M) 




Donor 


Classification 


A/T f 




[days] 


[M e ] 




Spect. Type 




l M el 


GRS 1915+105 a 


33.5 


9.5 


± 3 





K/M III 


LMXB/Transient 


14 ± 4 


V404 Cyg 


6.471 


6.09 


± 


04 


K0 IV 


1 1 


12 ± 2 


Cyg X-l 


5.600 


0.244 


± 


005 


09.7 lab 


HMXB/Persistent 


10 ± 3 


LMC X-l 


4.229 


0.14 


± 


05 


07 III 




> 4 


XTE J1819-254 


2.816 


3.13 


± 


13 


B9 III 


IMXB/Transient 


7.1 ± 0.3 


GRO J1655-40 


2.620 


2.73 


± 


09 


F3/5 IV 


J7 


6.3 ± 0.3 


BW Cir 6 


2.545 


5.74 


± 


29 


G5 IV 


LMXB/Transient 


> 7.8 


GX 339-4 


1.754 


5.8 


± 


5 


_ 






LMC X-3 


1.704 


2.3 


± 


3 


B3 V 


HMXB/Persistent 


7.6 ± 1.3 


XTE J1550-564 


1.542 


6.86 


± 


71 


G8/K8 IV 


LMXB/Transient 


9.6 ± 1.2 


4U 1543-475 


1.125 


u.zo 


it u 


01 


A2 V 


IMXB/Transient 


9.4 ± 1.0 


H1705-250 


0.520 


4.86 


± 


13 


K3/7 V 


LMXB/Transient 


6 ± 2 


GS 1124-684 


0.433 


3.01 


± 


15 


K3/5 V 


)? 


7.0 ± 0.6 


XTE J1859+226 c 


0.382 


7.4 


± 1 


1 




71 




GS2000+250 


0.345 


5.01 


± 


12 


K3/7 V 


?7 


7.5 ± 0.3 


A0620-003 


0.325 


2.72 


± 


06 


K4 V 


5) 


11 ± 2 


XTE J 1650-500 


0.321 


2.73 


± 


56 


K4 V 


77 




GRS 1009-45 


0.283 


3.17 


± 


12 


K7/M0 V 


77 


5.2 ± 0.6 


GRO J0422+32 


0.212 


1.19 


± 


02 


M2 V 


77 


4 ± 1 


XTE J 11 18+480 


0.171 


6.3 


± 


2 


K5/M0 V 


77 


6.8 ± 0.4 



f Masses compiled by |Orosz (2003) | and [Charles fc Coe (2006)| 

a New photometric period of 30.8+0.2 days recentfy reported by |Neil, Bailyn fc Cobb (2006) | 
The implied ma ss function, assuming constant velocity amplitude, would be 8.7 Mq. 
6 Updated after |Casares et al (2007)] 



Period is uncertain, with another possibility at 0.319 days (see Zurita ct al 2002). This would 
drop the mass function to 6.18 Mq. 



ordered by orbital period. The case of GRS 1915+105 is noteworthy, not only because 
of its long orbital period and large mass function; also because IR spectroscopy was 
needed to overcome the > 25 magnitudes of optical extinction and reveal the radial 



velocity curve of the companion star (Greiner, Cuby & McCaughrean 2001). It should 



be noted that, although some systems have mass functions < 3 Mq, solid constraints on 
the inclination and/or M c can be set which result in M x > 3 Mq. The great majority of 
BHs are transients and only 3 show persistent behaviour, the HMXBs Cyg X-l, LMC X-l 
and LMC X-3. From the list we also note that a new class of transient X-ray binaries with 
A-F companions is starting to emerge, the so-called Intermediate Mass X-ray Binaries 
(IMXBs) . It has been proposed that LMXBs may descend from IMXBs through a phase 
of thermal mass-transfer ( |Pfahl, Rappaport &: Pod siadlowski 2003), although it is still 
unclear whether IMXBs and LMXBs represent an evolutionary sequence. 

3.1. A novel technique: fluorescence emission from the irradiated donor 

In addition to the dynamical BHs, there are ~ 20 other BH candidates based on their 
X-ray temporal and spectral behaviour (McClintock & Rcmillard 2006). Unfortunately, 
they have never been seen in quiescence, or they simply become too faint for an optical 
detection of the companion star. However, a new strategy was devised to allow the 
extraction of dynamical information in these systems during their X-ray active states. It 
utilises narrow high-excitation emission lines powered by irradiation on the companion 
star, in particular the strong CIII and fluorescence NIII lines from the Bowen blend at 
AA4630-40 (see Fig. 3). This technique was first applied to the NS LMXB Sco X-l and the 
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Dopplcr shift of the CIII/NIII lines enabled the motion of the donor star to be traced 
for the first time ( |Steeghs fc Cas arcs 2002 ). This was also attempted during the 2002 
outburst of the BH candidate GX 339-4, using high-resolution spectroscopy to resolve 
the sharp NIII/CIII lines. The right panel in Figure 3 shows the radial velocity curve of 
the NIII lines folded on the 1.76 day orbital period. The orbital solution yields a velocity 
semi-amplitude of 317 km s _1 which defines a strict lower limit to the velocity amplitude 
of the companion star because these lines arise from the irradiated hemisphere and not 
the center of mass of the donor. Therefore, a solid lower bound to the mass function is 5.8 
M© which provides compelling evidence for a BH in GX 339-4 ( |Hynes et al 2003] ) . And 
most important, this technique opens an avenue to extract dynamical information from 
new XRTs in outburst and X-ray persistent LMXBs, which hopefully will help increase 
the number of BH discoveries. 




Figure 3. Detecting the companion star in X-ray active LMXBs. Left: the main high excitation 
emission lines due to by irradiation of the donor star in Sco X-l. Adapted from Steeghs & 
Casares (2002). Right: Radial velocities of the NIII lines (top), Hell A4686 (bottom) and the 
wings of Hell A4686 (middle) in GX 339-4 folded on the 1.76 day orbital period. After Hynes 
et al (2003). 



3.2. Further BH signatures 

Aside from dynamical arguments, there are other observations that lend support to the 
absence of a solid surface and, hence, the BH nature of these accreting compact objects, 
namely: 



• Lack of pulses and Type I X-ray bursts. This has been quantified by Remillard et al (2006) 



using observations of dynamical BHs over 9 years of RXTE data. The probabiliy that 
the non-detection of bursts were consistent with a solid surface is found to be ~ 2 x 10~ 7 . 

• The classic colour-colour diagram of X-ray binaries at high accretion rates shows a 
clear separation in the evolution of NS and BH binaries (jDone fc Gicrlihski 1999|. This 
has been ascribed to the presence of a boundary layer in NS which gives rise to an 
additional thermal component in the spectrum and drags NS outside the BH region. 

• For a given orbital period, quiescent BH binaries are ~ 100 times dimmer than quies- 
cent NS binaries (Me nou et al 1999a|) . This difference is interpreted as thermal radiation 
from the NS surface which, for BHs, is advected through the event horizon. 
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4. BH demography 

Despite the handful of confirmed BHs one can try to account for selection effects and 
estimate the size of the underlying galactic population. To start with, dynamical studies 
of XRTs indicate that about 75% contain BHs, i.e. M x > 3 Mq. Also, the extrapolation 
of the number of BH XRTs detected since 1975, with outburst duty cycles of ~ 10 — 100 
years, suggests that there is a dormant population of ~ 10 3 BH binaries (Romani 1998 
and references therein) . This is likely to be an underestimate if one accounts for systems 
with longer recurrence times and a likely population of faint persistent BH LMXBs 
(jMenou et al 1999b). Even so, incidentally, these numbers are in reasonable agreement 
with recent population synthesis calculations of BH binaries ( |Yungelson et al 2006[ ) . On 
the other hand, stellar evolution models predict a population of about 10 9 stellar-mass 
BHs in the Galaxy ( Brown & Bethc 1994). Therefore, our observed sample of BH XRTs 
is just the tip of the iceberg of a large hidden population of which nothing is known. Is 
it still possible to extract some meaningful information from such a limited sample? 




Figure 4. Left: Mass distribution of compact objects in X-ray binaries. Arrows indicate 
lower limits to BH masses. Right: observed mass distribution of compact objects in X-ray 
binaries (shaded histogram), compared to the theoretical distribution computed in Fryer 
fc Kalogera (2001) for the "Case C + Winds" scenario (dotted line). Mass- loss rates by 



Woosley, Langer & Weaver (1995) were used in the computations. The model distribution has 
been re-scaled for clarity. 

The main property of the BH population is the mass distribution. However, in order 
to get accurate BH masses, in addition to the mass function one needs to determine the 
binary inclination and the mass of the companion star or a related quantity such as the 
mass ratio. This information can be extracted from two experiments: (i) resolving the 
rotational broadening V TO t sin i of the companion's absorption lines, which is correlated 
with the binary mass ratio and (ii) fitting synthetic models to ellipsoidal lightcurves 
from which the inclination angle can be determined. This is the classic method to derive 



masses which has been reviewed in several papers e.g. Casares (2001) 



Following this prescription, reliable BH masses have been determined in 15 binaries. 
These are listed in the last column of Table 1 and displayed in Figure 4, relative to 



NS masses compiled by Stairs (2004) and Lattimer & Prakash (2004) The figure nicely 



shows a mass segregation between the 2 populations of collapsed objects, with NS clus- 
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tcring around 1.4 M Q and BH masses scattering between 4 and 14 M Q . Many more BH 
masses and errorbars ^10% are needed before fundamental questions can be addressed 
such as (i) do BHs masses cluster at a particular value? (ii) what are the edges of the 
distribution, i.e the minimum and maximum masses for BHs to be formed in X-ray bi- 
naries? (iii) is there a continuum distribution of masses between NS and BHs? All these 
questions are intimately related to models of SN and close binary evolution. 

Figure 4 also shows the histogram of the current distribution of compact remnant 
masses, relative to a model distribution computed by Fryer & Kalogera (2001) which 
includes mass-loss through winds and binarity effects. The model shows a mass cut at 
12 M© but this seems to be challenged by the compact objects in V404 Cyg and GRS 
1915+105, with 12 ± 2 Mq and 14 ± 4 Mq respectively. However, the conflict may 
be spurious since recent estimates of mass-loss rates in Wolf-Rayet stars suggests that 
previous determinations were biased too low ( Nugis k, Lamers 2000[ ). Also, the model 
predicts a continuum distribution between NS and BH whereas observations seem to 
show a paucity of objects at ~3-4 Mq. This may be caused by selection effects since low- 
mass BHs are expected to be persistent. Recently, compact object masses in the range 2-4 
M have been reported in 4U1700 -37 (jClark et al 2002|) . V395 Car (jShahbaz et al 2004|) 
and LS 5039 (jCasares et al 2005)) . Although they are assumed to contain NS, none of 
them has shown X-ray bursts nor pulses. Therefore, they may well be members of the 
missing low-mass BH population. We are clearly dominated by low number statistics and 
more observations are required. 



5. Conclusions 

The best observational evidence for the existence of stellar-mass BHs is provided 
by dynamical studies of X-ray binaries The first solid candidates were the classic X-ray 
binaries Cyg X-l and, in particular, A0620-00. BHs became confirmed with the discovery 
of mass functions in excess of 5-6 M , the first case being V404 Cyg. Their global X-ray 
properties (such as the lack of pulsations/bursts, weak quiescent Lx, etc.) also support 
the presence of an event horizon in these objects. 

XRTs are the best hunting ground for new stellar-mass BHs with 17 cases currently 
known and estimated masses between 4-14 M Q . These are are only the tip of the iceberg 
of an estimated dormant population of ~ 10 3 BH binaries and ~ 10 9 stellar-mass BHs 
in the Galaxy. Clearly many more discoveries and better statistics are essential to derive 
useful constraints on BH formation models. 

New strategies, aimed at unveiling more quiescent BH transients, and novel tech- 
niques, such as the detection of the irradiated donor using high-excitation reprocessed 
lines, need to be exploited. 
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Ramesh Narayan: When you make use of the ellipsoidal model to determine the binary 
inclination, how do you avoid contamination of the stellar light by the disc light? 
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JORGE Casares: The example I showed is a "textbook" lightcurve of J1655-40 where 
the hot F6 III companion star totally dominates the optical light. However, for the great 
majority of black hole binaries (which contain much cooler companions in shorter orbital 
periods) the contamination by the accretion disc light certainly starts to show up in 
the form of flaring activity and superhumps. It is commonly assumed that the active 
contribution of the accretion disc light decreases at longer wavelengths and, hence, IR 
lightcurves are used to estimate the orbital inclination. 

Tsevi Mazeh: For HMXB, is it a necessary assumption that the star fills its Roche 
lobe? If not, the ellipsoidal effect is difficult to model. 

JORGE Casares: Roche- lobe filling is not a general assumption of the models because 
HMXBs are normally wind-fed systems. Normally, the ellipsoidal fits are applied with an 
extra free parameter which is the star filling- factor. 

Xiaopei Pan: Physical parameters of black holes have huge uncertainties. The Space 
Interfcrometry Mission (SIM) can provide high-accuracy determination of inclination, 
mass, distance, etc. The precisions of black-hole masses can be directly determined to 
5%. The inclination precision can reach 2%, and the distance of black holes can be 
measured to better that 2% by the SIM project. We hope in more collaborations on 
black hole issues. 

JORGE Casares: The error budget in black- hole mass determination is clearly dominated 
by uncertainties in the inclination angle and indeed the improvement that the SIM project 
may provide will be extremely useful. 

Andrew King: Comment: I would like to reinforce your point about the black holes 
we detect being the top of an iceberg. Long-period LMXBs have wide separations, big 
discs, and are all transient. However, the interval between outbursts is so long that we 
will never detect them at all. We can detect the endpoints of these LMXBs when they 
contain neutron stars as very wide binaries with millisecond pulsars, but not when they 
contain black holes. 



